Optimal Phase Description of Chaotic OsciUators 
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We introduce an optimal piiase description of chaotic oscillations by generalizing the concept of 
isochrones. On chaotic attractors possessing a general phase description, we define the optimal 
isophases as Poincare surfaces showing return times as constant as possible. The dynamics of the 
resultant optimal phase is maximally decoupled of the amplitude dynamics, and provides a proper 
description of phase resetting of chaotic oscillations. The method is illustrated with the Rossler and 
Lorenz systems. 
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I. INTRODUCTION 

Phase description lies in the base of theory of self- 
sustained, autonomous oscillators [H-Ill. A prudently de- 
fined phase variable yields a one-dimensional description 
of the oscillator, allowing one to characterize important 
aspects of its dynamics such as regularity of oscillation, 
sensitivity to external forcing, etc. Moreover, the concept 
of phase is important for data analysis of oscillatory pro- 
cesses in physics, chemistry, biology, and technical appli- 
cations, where various approaches exist for extracting dif- 
ferent variants of phase variables from oscillatory scalar 
time series. 

On a very basic level, every phase description starts 
with the identification of those states of the oscillator 
which are in the same phase. For a good phase descrip- 
tion, the identification must be done in an invariant way 
- independent of the variables and observables used - 
in order to make statements about the oscillator's phase 
dynamics non-arbitrary and comparable. The standard 
procedure of phase reduction is valid for periodic oscil- 
lators that possess a stable limit cycle. There, a certain 
family of Poincare sections, called isochrones, is used for 
the identification of states: Each isochrone consists of 
those states which are mapped onto each other after one 
oscillation period T , and which converge to the corre- 
sponding state on the limit cycle [1, Q . 

Even though chaotic oscillators do not possess a sta- 
ble limit cycle, a phase-like variable has been used for 
their description. In this sense, the phase dynamics of 
chaotic systems has been initially discussed in relation 
to diffusion properties of phase 0, 0] and to phase syn- 
chronization fsl-UOll. However, to describe these features, 
one does not need a good microscopic, i. e. on the time 
scale of the order of a characteristic period T, definition 
of the phase because both diffusion and synchronization 
are defined macroscopically, i. e. for time scales much 
larger than T. On the other hand, in the theoretical de- 
scription of phase synchronization a proper microscopic 
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phase definition was pre-assumed [111 Il2l |. although no 
practical algorithm for construction of a phase variable 
with good properties has been presented. The reason is 
that the chaotic phase diffusion destroys a rigorous no- 
tion of time-coherent isochrones, because any two states 
of the chaotic oscillator which are thought to be initially 
at the same phase will diverge as their respective phases 
diffuse. 

In this article, we suggest a numerical technique for 
phase description of chaotic oscillations. For this goal, we 
construct special Poincare sections, which we call optimal 
isophases. The choice of isophases is based on the prop- 
erties of their return times. In an application to chaotic 
oscillators, we demonstrate an intimate relation between 
optimal isophases, chaotic phase diffusion, and unstable 
periodic orbits. Specifically, we discuss the reduced phase 
dynamics of chaotic oscillations and the decoupling of the 
amplitudes from the phase dynamics. Next, we use the 
optimal phase to introduce a proper framework for the 
description of phase resetting of chaotic oscillators. 

Starting with an outline of the standard phase defi- 
nition for periodic oscillators via the isochrones, we in- 
troduce the generalized concept of isophases of chaotic 
oscillators in Sec. HT] In Sec. IIIIl certain dynamical prop- 
erties of the optimal phase are highlighted by the ex- 
ample of the Rossler oscillator. Thereafter, the relation 
between optimal isophases and unstable periodic orbits 
is presented (Sec. IIVI) . In Sec. |Vl certain aspects of the 
theory are presented for the Lorenz oscillator. In the last 
section we discuss our results. 



II. ISOPHASES OF PERIODIC AND CHAOTIC 
OSCILLATORS 

A. Periodic oscillators and their isochrones 

Phase is a natural variable for the description of peri- 
odic motions in dynamical systems. It can be introduced 
in different ways, with different levels of mathematical 
rigor di, [1, 13, Here we outline an approach that is 
mostly suited for a generalization to the case of chaotic 
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systems. 

The consideration starts with a general dissipative 
dynamical system showing stable periodic oscillations; 
the system's state x(i) is, thus, attracted to a limit 
cycle xo(t) having period T. In a vicinity of this 
periodic attractor the state space can be foliated by 
a non- intersecting family of Poincare sections J(<p), 
parametrized by a phase variable ip with period 27r. With 
J{'f), a phase variable Lp{t) can be assigned to each state 
of the trajectory x(t) g J{(p{t)). Therefore, the fam- 
ily of isochrones J(<p) provides a precise definition of 
what is meant by an oscillation: The system completes 
one oscillation if the variable (p grows by 27r, i.e. if the 
trajectory returns to a chosen isochrone, consequently 
passing through all sections in J{(p). In order to sim- 
plify nomenclature and to distinguish this variable from 
the genuine phase, introduced below, we term ip as pro- 
tophase. Introducing coordinates on the sections Ji^p), 
one can parametrize each point by a vector of amplitudes 
a and the protophase ip. 

There are various equivalent ways to foliate the state 
space in such a way that ^p grows monotonically; for pe- 
riodic oscillators with a period T, the optimal foliation 
does exist [1]. It can be introduced by considering the 
stroboscopic map x(t) — > x(i -I- T). Clearly, all points 
on the limit cycle are stable fixed points of this map. 
Hence, for each fixed point Xq there exist a stable mani- 
fold which converges to xq under the action of the stro- 
boscopic map. These stable manifolds, called isochrones, 
constitute a special foliation of the neighborhood of the 
limit cycle, for which by construction the Poincare map 
is the same as the stroboscopic map. 

In this way one introduces the phase of oscillation so 
that its time evolution does not depend on the amplitudes 
a. By virtue of a trivial reparametrization ip ^ 9 — 

/ '^^ ^^^^ foliation, one can introduce the genuine 
phase 9 which grows strictly uniformly in time, with a 
constant instantaneous frequency 9 = u = 27r/T. This 
phase, defined in the whole basin of attraction of the 
limit cycle, serves as a basis for a theoretical description 
of perturbed periodic oscillations fll|. In particular, one 
can easily formulate phase-resetting properties in terms 
of this phase: If a state on the limit cycle x' is instantly 
perturbed to some other state (even outside of the limit 
cycle), x' — > x" , then the phase is reset by a value 
A9 = 6'(x") — 6'(x') , which remains constant in the course 
of further evolution (see also Sec. IIII D] below). 

Noteworthy, the extension of the phase to a vicinity of 
a periodic orbit can be defined either for a stable or un- 
stable limit cycle. In the latter case, instead of using the 
stable manifold, one constructs the isophases by using the 
unstable manifolds of the fixed points of the stroboscopic 
map. However, for saddle limit cycles having both sta- 
ble and unstable directions, this construction fails. Here 
one can construct isochrones on the stable and unstable 
manifolds separately, but not in the whole vicinity of the 
cycle. With this in mind, we use below for the chaotic 
case, where isochrones do not exist, the term "isophases" 



instead of the usual "isochrones" . 



B. Protophase for chaotic oscillators 

We start the generalization of phase description to 
chaotic oscillators by discussing the construction of the 
protophase. For this purpose, we need the chaotic at- 
tractor to show the same property as a limit cycle, 
namely that there exists a family of non-intersecting 
Poincare sections Ji^p), monotonically parametrized by 
a protophase p. The requirement includes periodicity, 
J{ip + 27r) = J{'~p), and that any trajectory on the at- 
tractor successively crosses each Poincare section J{p>) 
transversally. Of course, not all chaotic attractors pos- 
sess such a family, but those which have such a foliation 
can be described in terms of phases and are the subject 
of further consideration here. 

Let us consider as an example the Rossler oscillator (Tsj 

i = -y - z , 

ij^x + O.lby, (1) 
i = 0.2 + z(a;- 10) , 

and take a family of Poincare sections Ji^fi) defined via 
the cylindrical coordinates 

pi = tan^i - ; a = (r, h) = (^y/x^+y^, . (2) 

This family of Poincare sections with constant pro- 
tophase (fi is shown in Fig.[lja)). However, other families 
can be defined as well; an example of another foliation 
based on the protophase <^2 = Vi + 0.7 In r is counter- 
posed in Fig. mb). 

Because the difference of any two protophases is 
bounded, the asymptotic properties of their phase dy- 
namics, such as the mean frequency and the diffusion con- 
stant of the phase rotations, do not depend on the defini- 
tion of the protophase. However, local, microscopic prop- 
erties of the dynamics for two protophases are different 
as it becomes apparent through the irregularly fluctuat- 
ing phase difference ipi{t) — (p2{t) shown in Fig.[ljc). The 
fluctuations show a bounded and irregular pattern that is 
speciflc to arbitrarily chosen variants of the Poincare sec- 
tions. In order to define a "genuine" phase, such as that 
of periodic oscillators, we need to define the "isophases" 
of chaotic attractors. (Notice that for chaotic oscillators 
the isochrones generally do not exist.) Because the phase 
of a chaotic system is in fact not as "genuine" and unique 
as in the periodic case (see discussion below), we will re- 
fer to it as optimal phase, in the sense that it represents 
oscillating properties of chaos in an optimal way. 

C. Optimal isophases for chaotic oscillators 

The genuine phase of periodic oscillators is defined 
by the basic property that there exist Poincare sections 
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FIG. 1. (Color online) (a,b) Two different families of Poincare 
sections or the Rossler system (red bold lines). Both families 
yield a proper definition of oscillation. The corresponding 
protophases <^i,2 are, however, different, so that 'P2{t) — ipi(t) 
shows irregular bounded fluctuations (c), specific to the par- 
ticular shapes of the Poincare surfaces. 

where all return times are exactly equal to the period 
of oscillations; i. e. the corresponding Poincare maps are 
stroboscopic maps as well. Naturally, such a situation 
does not generally occur for chaotic oscillators. On the 
one hand, this is plausible because different periodic or- 
bits embedded in chaos usually have different basic peri- 
ods (total period divided by the number of crossings with 
a Poincare surface, see Eq. ([T]) below). On the other 
hand, a coincidence of Poincare and stroboscopic map 
would also imply the absence of phase diffusion what, 
however, is a degenerate, sparsely observed situation 

Since isophases of chaotic oscillators defined as sec- 
tions with constant return times do not exist in the strict 
sense, we introduce optimal isophases that approximate 
the property above with some accuracy. Practically, we 
construct the optimal isophases as a smooth Poincare sec- 
tion with a minimal (bounded by the smoothness) varia- 
tion of return times. As this condition is not unambigu- 
ous, we describe below an algorithm that we practically 
use. 

The starting point of our construction is a suitable 
vector time series x(t), < t < tend, of chaotic dynam- 
ics, which can be obtained by numerical simulation or 
by embedding observed oscillations [l5|. The first step 
is to introduce an arbitrary protophase ip as described 
above. Using it, we can estimate the average period of 
oscillations as 

^{tend) - </j(0) 

With this period, we define a family of stroboscopic sets 
for the trajectory x(t) as 

x,(0~) = x(^r + fcr) , k^i,2,...,K,nd- (3) 



Here 6 ^ [0, 2tt) serves as a (still preliminary) phase 
parametrizing stroboscopic sets, and each set consists of 
Kend points. Thcsc sets are invariant under the stro- 
boscopic map with time interval T, but cannot serve as 
Poincare maps as they are not smooth curves, because 
the rotation in chaotic systems is non-uniform. The 
larger the total time interval t^nd, the stronger is the 
spreading of the points of the stroboscopic set. We illus- 
trate this in Fig. [2j We note, that only in a degenerate 
case where the phase diffusion of the chaotic oscillator 
vanishes, these stroboscopic sets would be smooth lines 
which can be used as Poincare sections; such degenerate 
chaotic attractors (see an example in [l^) possess the 
same rigorous phase description as periodic oscillators. 
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FIG. 2. (Color online) The stroboscopic sets [Eq. ((3|] for 
the Rossler attractor for two lengths of the trajectory (a) 
tend = 10'' and (b) tg„d = 5 • 10'' are shown by squares. The 
trajectories are shown with gray lines, the optimal isophases 
obtained by fitting the by a polynomial tf(r) of order four are 
shown with black bold lines. 

In order to obtain a proper smooth Poincare section, 
we fit the stroboscopic set, in the sense of least squares, 
by a polynomial ip = i p(a') (we use standard fitting pro- 
cedure as described in [l7[). The resulting curves shown 
in Fig. [5] are our optimal isophases, i.e. the curves of 
constant phase 9. 

If we restrict ourself to rather smooth isophases only, 
a good practical approximation can be achieved if one 
introduces a global phase correction function A according 
to 

= ^ + A(^,a) (4) 

and finds its representation in terms of polynomial basis 
functions: For each of the amplitude components aj we 
use the powers a", and for the phase variable ip we use 
trigonometric polynomials exp(i(/?Z). For example, for the 
Rossler system in 1-1-2 dimensions, consisting of phase 
radius r and height ft,, the phase correction is represented 
using a set of cocfhcicnts CmnV- 

Nr Nh 

A{ip, r,h)=Y,Y.Yl r"ft"e''^ . (5) 

m=0 n=0 1=0 

The coefficients can be computed by applying a linear 
least squares fit [l3| to the stroboscopic sets. In this 
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FIG. 3. (Color online) A global approximation of optimal 
isophases (blue dots which look like bold lines), obtained for 
the Rossler attractor (gray) using the approximation ((5| with 
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FIG. 4. (Color online) Return times r„ for the Rossler oscilla- 
tor [Eq. IT])]. Filled black squares correspond to an arbitrary 
Poincare section y — 0,x < 0, here the spreading of the re- 
turn times is large. Local (blue open circles) and global (red 
crosses) approximations (nearly coinciding on the figure) of 
the optimal isophases yield a strongly reduced spreading of 
the return times. 



way it is easy to find an optimal phase globally, as a 
function of the state space coordinates x. We illustrate 
the isophases obtained in this way in Fig. [3l 

In Fig. |4] we compare the quality of the optimal 
isophases obtained via representation ([S]) with the results 
of the local fitting of stroboscopic sets as in Fig. [2] We 
compare the return times for these isophases with the re- 
turn times of the Poincare section y = 0,x < 0. One can 
see that globally defined smooth isophases in the form 
([5]) give a quite good minimization of the variability of 
return times. 



III. DYNAMICS OF THE OPTIMAL PHASE 

In this section we discuss dynamical properties of the 
optimal phase introduced with help of optimal isophases. 



A. Return time map 

A natural way to characterize the time intervals T„ 
between successive crossing of a Poincare surface is to 



FIG. 5. The return time map [Eq. ^] of the Rossler oscil- 
lator [Eq. ([1])] for the cylindrical Poincare section ip = iir/S 
[Eq. ((2])] can be described as a one-dimensional chaotic map 
(black dots in (a)). Using the optimal isophase one obtains a 
map in a much smaller range (small box in (a) is enlarged in 
(b)). 



construct the return time map 

T„+i = M(T„) 



(6) 



In fact, because r„ is a function of the Poincare map 
coordinate, it is just a scalar observable, and M(r„) is 
not a function but rather a one-dimensional projection of 
a Cantor set. Nevertheless, for nearly two-dimensional 
strange attractors the Poincare map is nearly one- 
dimensional, and ^ looks like a curve (see Fig. [5K)- In 
Fig. Owe demonstrate, how this return time map changes 
if one uses an optimal isophase as a Poincare surface. 
First, the range of variations of r„ drastically shrinks. 
Second, one can hardly recognize the one-dimensional 
structure of the map: because now T„ is a "bad" observ- 
able, it does not reproduce the nearly one-dimensional 
nature of the Poincare map a„ — >■ a„_|_i. This means, 
that with the dynamics of the new optimal phase looks 
like a random process even on a microscopic time scale, 
that is of the order of the period T. 



B. Uniformity of phase rotations 

The basic property of the phase for a periodic oscilla- 
tor is that it rotates uniformly. For the optimal phase 
of a chaotic oscillator we cannot expect pure uniformity, 
but nevertheless it should be considerably increased com- 
pared to an arbitrary protophase. We illustrate this in 
Fig. \6\ Here we show the velocities of the protophase 
ip defined according to Eq. ((2]) and that of the optimal 
phase 9 defined according to isophases shown in Fig. |31 
While fluctuations in the protophase velocity ip heavily 
depend on ip, the fluctuations of 9 are almost uniformly 
distributed and, notably, in some regions are larger than 
those of the protophase. Similar results are reported in 
Ref. 
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We conclude that optimal isophases not only 
eliminate the amplitude dependence of the phase velocity, 
but they also flatten the phase dependence of its velocity 
fluctuations. 
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FIG. 6. The phase velocities for the Rossler system for the 
protophase defined according to ((2| (panel (a)) and according 
to isophases Fig. [31 panel (b). 



C. Decoupling of amplitude and phase dynamics 

One of the goals of introducing the phase is to decou- 
ple its dynamics from that of the amplitude. For periodic 
oscillators this decoupling is perfect, whereas for chaotic 
oscillators, it is only approximate. To illustrate, how cor- 
relations of the phase dynamics with the amplitudes are 
reduced when the optimal phase is introduced, we per- 
formed a "mixing" experiment; the results are depicted 
in Fig. [71 We started an ensemble of initial conditions 
on a certain Poincare surface and followed them for a 
time interval of length 5T (five average rotation peri- 
ods). The trajectories starting at small, medium, and 
large amplitudes a are marked separately in the plot. In 
the upper panels (a,b), where a "cylindrical" Poincare 
surface 0i = const is used, we see that the states started 
at small amplitudes are lagged behind while those started 
at large amplitudes are advanced. Contrary to this, us- 
ing the optimal isophase as an initial condition, we see 
that after 5 rotations all points are mixed and one can 
hardly distinguish the points which had different ampli- 
tudes at the beginning. This is another illustration of the 
fact that the dynamics of the optimal phase is effectively 
decoupled from the amplitude. 



D. Phase resetting of chaotic oscillators 

A basic application of the phase description of peri- 
odic oscillators is quantification of the system response 
to pulse stimulation by means of phase response curves. 
Given a state on the limit cycle x(0) one can determine 
the phase shift due to change of the state x{9) — > x' = 
x(6') -I- k simply by calculating 9' = 6'(x'). Because the 
phase rotates uniformly also outside of the limit cycle, the 
phase shift 6' — 9 remains invariant and characterizes the 
phase resetting (for noise-induced oscillations this notion 



FIG. 7. (Color online) Two "mixing" setups where initial 
states of the Rossler oscillator [Eq. H])] (marked as symbols) 
are chosen either on the Poincare section 0i = 47r/3 [Eq. ([2|] 
(a) or on the optimal isophase (c). Panels (b) and (d) show the 
same points at time t = ^ - T, where T is the average oscilla- 
tion period. States on the optimal isophase show less diffusive 
broadening in direction of the phase than the points on ar- 
bitrary Poincare section. Moreover, states of different ampli- 
tude become indistinguishable only for the optimal isophase, 
as seen by the mixing of markers. 



can be also introduced in an optimal sense [l9|). 

This approach has to be slightly modified when ap- 
plied to chaotic oscillators. If both states x and x + k 
lie on the attractor, then their optimal phases are well- 
defined, and the phase shift can be simply calculated 
as 9{-x. + k) — ^^(x). However, generally the state x + k 
lies outside of the attractor, and we have to generalize 
the definition of the optimal phases from the attractor 
to its vicinity. This is ambiguous, because the opti- 
mal isophases are not genuine isophases. They are not 
strictly invariant under time shifts and we cannot define 
the phase of the state x -|- k by following its time evolu- 
tion for arbitrarily large times. Instead, we have to fix 
the time interval after which the phase of the state x + k 
is defined. For the Rossler model we choose the mean 
period T as such an interval, as the relaxation time of 
approaching the attractor is typically smaller. So, we 
define 6'(x -I- k) = 9\T{-s. + k)] where T is the operator 
of time evolution over the average period T (see Fig. [5]). 
Applying now the representation ([5]), we obtain the PRC 
of the Rossler attractor PRC(k, x) 6'[f (x + k)] - 6'(x) 
as shown in Fig. [H 



IV. ISOPHASES AND UNSTABLE PERIODIC 
ORBITS 

In this section we discuss a relation between optimal 
isophases of a chaotic system and unstable periodic orbits 
(UPOs) Xo(i + t) = Xo(t) embedded in chaos. For each 
UPO one can define the phase on this orbit just from the 
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FIG. 8. (Color online) A state x of the Rossler attrac- 
tor (gray trajectory) is kicked to the state x + k (black 
points with an arrow). After one period the perturbed tra- 
jectory (dashed line) returns to the attractor and now lies 
on the isophase ^[^(x-l-k)] (red dotes). The kick's effect 
on the oscillator's phase is therefore given by the phase shift 
PRC(k, x) = 6l[f (x + k)] - 6»(x). 
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FIG. 9. (Color onhne) The phase resetting PRC(k, x) for the 
states on the Rossler attractor, color (gray scale) coded, for 
k= (1,0,0). 



condition of uniform rotation. This approach is discussed 
in Sec. IIV Al Similar to the construction discussed in 
Sec. Ill Al we can extend the notion of the phase for each 
periodic orbit to its stable or unstable manifold; these 
ideas are presented in Sec. IIVBI 



A. Approximation of orbit phase sets 

For each UPO one can introduce a topological period 
(lap number) p as the number of intersections with a 
Poincare section. With this number p and the total pe- 
riod r, we define the oscillation period 



z/ p 



(7) 



which is expected to be close, but not identical, to the 
mean period of chaotic oscillations (mean return time of 
the Poincare map). Next, for the UPO we can introduce 
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FIG. 10. (Color online) (a): The optimal isophase (black 
line, obtained using a fit with a high-degree polynomial) of 
the Rossler oscillator [Eq. ([1])] overlapped with orbit phase 
sets [Eq. ((§])] for a 9-orbit (green circles) and a 10-orbit (red 
squares), (b): A distance measure d [Eq. ([9])] quantifies how 
close is the orbit phase set to the isophase, here shown for 80 
p-orbits with p < 10. One can see weak correlations to the 
instability of the UPOs measured by the Floquet multiplier 
\P\- 



the phase 6 that rotates uniformly with frequency 2Tr/iy 
so that 9{t) = 6{0) + 2Tip. With help of this phase, a 
family of point sets I {9), called orbit phase sets, can be 
defined as points which are attained at constant time 
intervals, equal to the oscillation period S: 



I(e) = {^a{nS) I n = 0,...,p-l} 



(8) 



with some arbitrary choice of the zero phase. 

Let us now take a Poincare surface which passes 
through the orbit phase set I {9). (Of course, there are 
many possibilities to draw such a surface, e.g., one can 
use splines.) Then it will be an approximation to an 
optimal isophase, as at least on the orbit phase set all 
the return times will be equal to S. We illustrate this 
with Fig. [Tola), where we show orbit phase sets of two 
UPOs, with topological periods p = 10 and p = 9, for 
the Rossler system ((TJ . Since the orbits do not share any 
state, the zero phases can be chosen separately. Practi- 
cally, the phase offsets have been chosen in a way that 
the orbit phase sets are mostly close to each other and 
approximate the same isophase, which is also drawn for 
comparison. One can see that the orbit phase states in- 
deed can serve as approximations for the isophases. 

This approximation is expected to work better for 
larger periods and for periodic orbits which are most 
"typical" . The probability for a trajectory to approach 
the orbit depends on the stability of the UPO, quanti- 
fied by its unstable Floquet multiplier [2^. Therefore, 
it is expected that the correspondence between isophases 
and the orbit phase sets will be better for UPOs which 
are visited more often because they are less unstable. To 
check this for the Rossler system, we introduce a measure 
d of the distance of a p-orbit y to the optimal isophase 
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not close to the optimal isophases obtained by another 
methods. 

It is instructive to try to construct the isophase on the 
unstable manifold of the UPO using not its period S, but 
the mean period T. It is clear that such an isophase can- 
not exist, but trying approximate it (see appendix for 
details) we obtain a singular curve (Fig. [TT]) . This is an- 
other representation of non-smoothness of stroboscopic 
sets that appears in the algorithm described in Sec. Ill CI 
due to non-existence of true isophases. In fact, when one 
tries to construct an isophase, such a singularity will ap- 
pear for every periodic orbit, and the procedure should be 
constrained by the requirement that the isophase should 
be sufficiently smooth. 



FIG. 11. (Color online) For a UPO of the Rossler oscillator 
(dashed line), the orbit phase set [Eq. ^] can be extended 
to the unstable manifold in two ways. Blue solid line shows 
the extension where the orbit's period S is used. Red circles 
depict the extension based on the mean period of the chaotic 
attractor. The signularity of the latter curve indicates a di- 
vergence of the phase correction (see Appendix [A] for the an- 
alytic form of this correction for the unstable Stuart-Landau 
oscillator.) 



shown in Fig. |4] with a green line, as 



p-i 

'-El 

k=0 



(9) 



where are coordinates of the orbit phase set (for which 
we also optimized the zero phase to achieve a minimum 
of d) and yfe are the crossings of the periodic orbits with 
the isophase. This measure was calculated for the 80 
available UPOs together with their Floquet multipliers. 
It was found that orbits showing a larger distance had a 
tendency to be less stable (cf. Fig. ITUTb')'). 



B. Orbit isophase 



As described in Sec. Ill Al after the phase on a periodic 
orbit is introduced, the isophases in its vicinity can be 
defined separately on the stable and unstable manifolds 
of the orbit, as the stable and the unstable manifolds of 
the fixed points of the stroboscopic (with the period of 
the orbit) map. This definition can be applied to the 
UPOs in chaos, where the unstable manifolds are espe- 
cially interesting as they lie in the attractor. 

Let us consider the simplest UPO of the Rossler os- 
cillator that has topological period p — 1. Its oscilla- 
tion period is S" ~ 6.024 whereas the mean period of a 
typical trajectory is T « 6.073. Numerically, we calcu- 
lated the isophase on the unstable manifold of this orbit 
using the oscillation period S for the stroboscopic map 
and obtained the blue line in Fig. [TT] This isophase be- 
comes folded together with the unstable manifold, and is 



V. PHASE OF THE LORENZ SYSTEM 

In our presentation above we have used the Rossler 
model [Eq. ([!])] as the basic example. Here we discuss 
how the approach works for the Lorenz system 



X = 10 ■ {y — x) , 
y = 28 ■ X ~ y — xz 
8 

z — ~ — ■ z + xy . 



(10) 



Chaotic phase diffusion of the Lorenz system is orders 
of magnitude stronger than that of the Rossler oscillator 
Eq. ([1]), thus introducing its phase is a more challenging 
task. The main difficulty lies in the unboundedness of 
the return times of the Poincare map, due to presence of 
the saddle steady state at the origin {x — y = z = 0). 
Due to this, the stroboscopic sets are spread over the at- 
tractor and cannot serve as a basis for construction of 
isophases like described above. Therefore we applied the 
following iterative procedure for obtaining smooth opti- 
mal isophases. First, we use projections of the trajectory 
onto the plane (m = a/x^ + y^, z). On this plane the tra- 
jectory rotates around a center approximately at (12, 27), 
and the protophases can be easily defined (cf. [2l|). We 
choose a Poincare surface and find the points of the tra- 
jectory at the intersection with this surface, these are 

x{tk),y{tk), z{tk), k = 1,2, Of course, the times tk 

are not equidistant because the Poincare map is far from 
the stroboscopic one. We adjust the times tf, trying to 
make them equal, by introducing a parameter s on which 
these times depend, and letting them evolve according to 



dtk 
ds 



dV{h,t2...) 
dtk 



V 



2 ^ 

k 



it 



k+1 



tk-Tf 



(11) 

where T is the average period. One can easily see that 
the "evolution" of tk according to Eq. (ITTj) leads to equal- 
ization of the intervals tk+i — tfe because of minimiza- 
tion of the Lyapunov function V . However, we "evolve" 
the times tk only for a finite interval of s, and obtain 
new times tk = tk{s). The new points x{tk),y{tk),z{tk) 
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form a new, distorted and singular Poincare section. We 
smoothen this set by applying a kernel technique |23] and 
obtain a smooth new Poincare section with more equidis- 
tant time intervals. We make several iterations of this 
procedure, and finally obtain the approximate smooth 
isophases as depicted in Fig. [I2l 



40 



^ 20 




^ ' ' ' ^ 

-40 -20 20 40 

x + y 

FIG. 12. (Color online) Optimal isophases (different mark- 
ers/colors depict different isophases) of the Lorenz attractor 
[Eq. (Uni)] (grey line). 

To characterize the quality of the introduced isophases 
for the Lorenz system, we plot the return times for an 
initial arbitrary Poincare section and for the obtained 
isophase in Fig. [T31 We see that the variations of the 
return times decrease only slightly, and the singularity 
(corresponding to the stable manifold of the origin) re- 
mains. 




0.4 ^ ' ' ' 1 

5 10 15 20 

z 

FIG. 13. (Color online) Return times for the Poincare section 
II = 12, z < 27 (filled circles) and for the optimal isophases 
resulting from its iterations (shown in Fig. [T2]with black filled 
squares) (open circles). The variations of the return times 
only slightly decrease. 

In Fig.[T3]we use orbit phase sets of UPOs of the Lorenz 
system to approximate isophases. Nine periodic orbits of 
the Lorenz system with topological period 6 are shown 
with grey line. By manually adjusting phase shifts of 
these orbits, it is possible to arrange the isophase sets 
(different markers) to build a set close to a curve (drawn 
manually as a black line) which can serve as an optimal 
isophase. The form of this curve is close to one of the 
optimal isophases presented in Fig. [T^l 




^ ' ' ' ^ 
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FIG. 14. (Color online) Building an isophase using nine 
UPOs of the Lorenz system with p = 6 (see text for details) . 

VI. CONCLUSION 

In summary, we have proposed a method of phase re- 
duction of chaotic oscillators by generalizing the concept 
of standard isophases (isochrones) . In the absence of a 
stable limit cycle, the definition of optimal isophases of 
chaotic oscillations is solely based on their return times. 
Because of non- vanishing chaotic diffusion and embedded 
unstable periodic orbits with different periods, isophases 
could only be obtained in an optimal, approximate way 
constrained by certain smoothness conditions. In the 
case of the Rossler attractor, where the phase diffusion 
is relatively small, we obtain the optimal isophases by 
smoothing the stroboscopic sets of a chaotic trajectory. 
For the Lorenz attractor, where phase diffusion is large, 
an iterative numerical scheme was proposed. Using the 
Rossler oscillator as an example, we have presented dif- 
ferent aspects of the phase dynamics. Specifically, the 
decoupling of the phase dynamics from the amplitudes, 
as well as a way to describe phase resetting of chaotic 
oscillators have been outlined. 

The theory of optimal isophases can possibly provide 
a refined understanding of emergent behavior of weakly 
coupled oscillating systems. For example, a theoretical 
phase description of weakly coupled limit cycle oscillators 
can be extended to ones of greater complexity, such as 
stochastic or chaotic oscillators (cf. [Ullllj). In this way, 
one can treat more realistic models of natural systems. 
Furthermore, the theory can easily be utilized in the anal- 
ysis of observed chaotic oscillations where the numerical 
scheme described above can be used to refine a prelimi- 
nary phase description. This can help to reduce certain 
systematic errors which may be present in phase-related 
quantities such as coupling strengths. 

The theory is easily adaptable for the analysis of non- 
linear oscillations with a random component (for theoret- 
ical approaches see, e. g. [23l - |25j ). Here, the return times 
to optimal isophases have to be interpreted in an aver- 
age sense. The corresponding results will be presented 



9 



elsewhere. 



solutions of the equation 
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Appendix A: Isophases of unstable Stuart-Landau 
oscillator 



To give an analytically tractable example of isophases 
of UPOs on their unstable manifold, let us consider the 
unstable Stuart-Landau oscillator governed by 

r — r(r^ ~ 1); — a ~ nr^ . (Al) 

It is exactly solvable: For the initial conditions r(0) = R 
and 1^3(0) = $, it has the well-known solution 



r{t) = 



'1/2 



(A2) 



ip{t) = (a - K)t - K\nr{t) + $ + kIhR . (A3) 

Oscillator (jAll) shows an unstable periodic orbit (UPO) 
with frequency uj — a ~ n. Depending on the initial 
conditions, it either performs decaying oscillations (for 
i? < 1), or diverges in finite time (for R > 1). 

As the characteristic period we first choose that of the 
UPO: S — In order to obtain a phase which rotates 
independently of r, we set 9 — ujt + ^ + kIh R. Inserting 
9 into Eq. (|A3[) . we find that optimal isophases I (9) are 



= if + Klnr 



(A4) 



For each ($, R) e I{9) the return time for 9 ^ 9 + 2n is 
equal to S because 9 — lo. This is the standard definition 
of the isophases. 

Alternatively, one may think of the unstable Stuart- 
Landau oscillator as of a rarely visited part of the state 
space of a bigger chaotic system which has a different 
characteristic frequency ^ = cjq = <^ + Aw. It means, 
average period T is different from the period S of the 
UPO. Therefore, the condition that the return time for a 
Poincare surface is equal to T cannot be fulfilled on the 
orbit. To fulfill the condition for states off the periodic 
orbit, we now seek a phase with the dynamics 9 = ujq. 
Therefore, we rewrite Eq. (|A3p in terms of uj^t: 



ip{t) = ujQt + $ + K In i? - K In r - Awt(r) 



(A5) 



Here, we need to rewrite time as a function of radius, 
using ((X2|) . We get 



t{r) = Un\r^ 



11 - Ini 



In- 



R 



(A6) 



Vl-i?2 

After the substitution, a uniformly rotating phase is 
given by 9 = LOot + <^ + n\nR + l^u]\n{^/l - R^/R). Com- 
paring the result with Eq. (j?]), we obtain the phase cor- 
rection as 



A(r) 



Aw) In r In Ir^ 

' 2 ' 



II 



(A7) 



While the return time is equal to T off the periodic orbit, 
the phase correction diverges as ln|l — r| in the limit 
r — )■ 1. Thus, the "isophase" is singular and winds itself 
infinitely often around the limit cycle. 
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